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1. Introduction 




In the present paper an attitude servo is defined as a control system 
which adjusts the attitude and angular velocity of the spacecraft so as to 
follow the reference which may be varying in time. Such systems are multi- 
dimensional multivariable, nonlinear, and nonautonomous with infinite sets 
of forcing functions and initial conditions . A procedure for designing 
such systems for spacecraft which are controlled by means of an arbitrary 
angular momentum exchange and storage device is presented in the paper. 

The discussion is restricted neither to small angle motion nor to single 
axis motion of the spacecraft. Indeed, the nonlinear effects arising from 
kinematics and gyroscopic action, the presence of bounds on angular momen- 
tum exchange rate and storage capacity of the controlling device, and vari- 
ations in system parameters are taken into account. The proposed control 
scheme forces the spacecraft into the desired attitude by generating torque 
about all three axes simultaneously. 

2. A Model of Attitude Servos 

Consider three right-handed orthonormal triplets of vectors s, a, and 
d whose common origin is the fixed point of rotation of the spacecraft. 

Let s and a be fixed in inertial space and spacecraft, respectiv^ely, and 
let d represent the reference. The servo input, output, and error will 

A 

be defined by direction cosine matrices: the input A^g defines d rela- 

tive to s; the output Ag^g 


defines a relative to s; the error 
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R = (superscript t denotes transpose) defines a relative to d. 

Servo null occurs when R = I, the identity matrix. 

The model shown in figure 1 is proposed for the description of an 
attitude servo in which the controller is an arbitrary angular momentum 
exchange and storage device. The body coordinates of the angular momentum 
exchange rate are assumed to be the controlling variables of the servo. 
External torq.ues are assumed to be negligible. 

STATE a * (Ads. 
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Fig. 1 A model of attitude servos. 

The states of the servo are imbedded in a 24- dimensional space whose 
points a are represented mnemonically by the Quadruplet whose elements 
are the input matrix, error matrix, body coordinates of angular velocity of 
spacecraft, and inertial coordinates of total angular momentum of the sys- 
tem. The subset 9 is the region of operation of the servo on which ortho- 
gonality and saturation constraints are satisfied. The latter arise because 
any practical controlling device has limited angular momentum storage capac- 
ity which will be assumed to be spherically bounded by 

'^smax + <3max”max = ^max> -3max maximum principal moment of 

inertia of the main body. 

The state equations are defined on 0. The first two are the kinematic 
equations of input and error, respectively (for any column x, S(x) is 
defined by the skew symmetric matrix whose entries in the upper triangle 
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are X3, and x^). Kinematic eq,uations automatically preserve the 

orthogonality of and R. The second pair of state eq.uations are the 

dynamic eq.uations. The controlling variable is Zg^ = -(d/dt) (Ag^ghg - J^^a) 
where the matrix represents the moment of inertia of the main body in 

body coordinates. The angular acceleration depends on the control law £ 
which is assumed in the form of a state feedback z(a), and on the pertur- 
bation function n which is assumed to depend on a and the perturbation 
variable A. The last dynamic equation is a consequence of the assumption 
that external torques are negligible^ and automatically preserves the norm 
of 

There are two forcing functions. One is the angular velocity of the 
reference whose d-coordinates are denoted by the matrix w^. The other is 
the perturbation variable A. They are restricted to the sets of vector 
functions of time which are spherically bounded by 
respectively. 

The set Z of admissible control laws and perturbation functions is 
defined by two conditions. One requiring the control to be spherically 
bounded everywhere on 9 accounts for the fact that a practical control- 
ling device has limited angular momentum exchange rate. The other requir- 
ing the time rate of change of the magnitude of angular velocity to be not 
greater than zero on the velocity boundary of 9 forces the servo to 
remain in 0. 


3. Global Description of Attitude Servos 

Suppose that a definite admissible control law and perturbation 
function are given. The question whether the servo is fast enough for the 
mission requirements of the spacecraft may be resolved as follows. 

It is a consequence of Euler’s theorem on rotations that an angle 
(error angle cp) and an axis (error axis c) determine R. Thus^ 

R = exp[cpS(c)]. Conversely, cp = arc cos[ 0 . 5 (tr(R) - 1)], and for 
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cp € (o, 7t), c = 0.5 CSC Cp(r23 - 1*32, Tsi - ri 3 , ri 2 - 3 ^ 21 )^* The error 
angle cp is a ina thematically convenient and intuitively appealing scalar 
representation of attitude servo error, and forms a basis of the notion 
of a response envelope introduced next. 

For fixed admissible control lav and perturbation function the state 
equations determine the history 2 . “ 't) ^ t > 0 of the error 

angle for each admissible initial state a, reference velocity and 

perturbation variable A. The response envelope 2 ** defined as the 
envelope of all such histories. That is, 2 ** = 9**(*t), t > 0, where 


cp**(t) 


max 
a € 0 


Liid 


max 

€ W, e A 


cp(t, 


a, A) 


fixed t 


An approximation, which is in the spirit of the Liapunov theory, of 
the response envelope may be obtained as follows. Consider the surface 
V(a) = v(t). The fimction V(a) is assumed to be given explicitly, to have 
a gradient with respect to 00 everywhere on 0, and to include 0 for 
some finite v(0). Let the expansion rate of the surface be defined by 
V = W(v), where 


■ {a € e:V(a) = v, !!w^!! < 


1} (VaV(a.)d) 


•*amax> h'mi 

An upper estimate 2u response envelope 2 ** that 

cp^(t) > cp**(t) for every t > 0, may be obtained by maximizing at each t 

the error angle cp on this moving surface thus. 


/. \ max / V 

= {a € e;V(a) = v(t)} 


4. A Set of Control Laws 


Euler’s theorem gives rise to the concepts of error angle and error 
axis. The first was used above to describe an attitude servo. The second 
is used next for control. The set of control laws having the form 
indicated by the following equation (l) is proposed 
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z(a) = -Ja[gi(cp)g2(l|walp)c + ga(cp)wa] 


(I) 


Qualitatively, a control of this form may be thought as a generalization of 
the familiar position plus rate feedback used in one-dimensional servos. 

The first term in the brackets in equation (l) acts along the error axis; 
the second term acts along the velocity vector. The scalar functions gi, 
g2, and g3 are restricted only by the two conditions defining Z. 

An upper estimate of the response envelope corresponding to any con- 
trol law of this form may be computed by means of the V-f unction given by 
the following equation. 


V(a) = J 


gi(x) + ng3(x)sin 


i]-ur 


2 

dx 


+ ^ sin £ c^Wa 


g2(x) ' 2 

(I is an adjustable parameter. The actual computation of an upper estimate 
may be performed in a three-dimensional space whose points are related to 
states in 6 by the function 


q(a) 



sin 


2 . 

2 



5 . Example 


Consider the Orbiting Astronomical Observatory (OAO). It is con- 
trolled by means of three identical orthogonally placed reaction wheels. 

The ith con5)onent of Zq^ is the negative of the ith motor shaft torque. 

is the moment of inertia of the spacecraft with locked wheels minus 
the moment of inertia of the wheels about their spin axes. In what follows, 
\iax ^max denote the bounds on angular momentum storage capacity 

of each wheel and torque capacity of each motor, respectively. 

Let the spacecraft be controlled by (2), below, which is a special 
case of (1). 


z(a) = -Jg 


^max 


2Jmax L 


sat(cp, cp )c + 




max 




( 2 ) 
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sat(cp, cpg) equals cp/cpg for cp < cpg and 1 for cp > cpg. The dimensionless 
constant cpg = VxW)- 

Upper estimates of several response envelopes were computed on a 
digital computer using the following V-function. 

V(a) . ^ (1 - cos l) * i <p, + slo I 

The following plots show the results in dimensionless form. Consequently, 
they apply to any spacecraft, with cpg = 0.1, which is controlled hy (2). 




(a) Nominal servo. (b) Spherical error in acceleration. 

Fig. 2 Spacecraft response. 

Figure 2a shows upper estimates of response envelopes for the nominal 
servo for various bounds angular velocity of the input. The 

case with b = 0 corresponds to the state regulator. The plots indicate 
that regardless of in which admissible state the servo is initially, the 
attitude error will be less than 6^ after about 5 units of time (28 minutes 
for the OAO) even when the input varies arbitrarily with one quarter of 
maximum angular velocity allowed for the spacecraft. 

Figure 2b is relevant to the design of angular momentum exchange and 
storage devices. For example, suppose that the spacecraft is to be con- 
trolled not by reaction wheels but by a set of control moment gyros in some 
configuration. If the angular momentum of each gyro is approximated by its 
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spin momentum, the angular momentum stored in the complete device may he 
expressed as a function h(x) of the ginibal angles x of the gyros. The 
control Zg, = -(Sh/Sx)x. Let the exchange law of the device he 
X = F(x)z(a). If z(a) is given hy ( 2 ), cpg = 0 . 1 , and the following 
inequality holds for all possible x, then the response envelope of the 
spacecraft is hounded from above hy the ciirve h = O-3. 


^max 

Jmin 


I 



= h < 0.3 


Two reasons why h may he greater than zero are desired simplicity of 
F(x) and failure of some gyros in the device. 




(a) Spherical error in error axis. (b) Gyroscopic effects. 


Fig- 3 Spacecraft response. 

Part (a) of figure 3 is relevant to the design of attitude sensors. 

For example, let the attitude of the spacecraft he measured with a set of 
star trackers. The gimhal angles x may be expressed as a function 
g(Ag^g, Xq) of the attitude of the spacecraft Ag^g and inertial coordinates 
Xq of the guide stars. Let the output y of the gimhal processor he some 
function f (x, Suppose that y is used in place of sat(cp, cpg)c 

in the control law ( 2 ). Then for cpg = 0.1 the response envelope of the 
spacecraft is hounded from above hy the curve h = 0.4 if the following 
inequality holds for all orthogonal R and all expected inputs and guide 


stars. 
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f(g(RAets, xj, Aas, Xq) - sat(cp, cpg)c 
sat(cp, cpg) 

Two reasons why b may be greater than zero are the desired simplicity 
of the processor f and failure of some trackers. 

Part (b) of figure 3 shows the influence of gyroscopic coupling 
arising from nonzero total angular momentum. The perturbation function 
n(a. A) defined in the figure generates a set of vectors which always 
includes the gyroscopic terms. It is seen that in the case of the OAO 
gyroscopic coupling is not very significant even when the system is loaded 
with as much as O.3 of its angular momentum storage capacity. 


= b < 0.4 



